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Abstract 



In this paper we revisit Brill's proof of positive mass for three-dimensional, 
' time-symmetric, axisymmetric initial data and generalise his argument in var- 

I ions directions. In 3 -I- 1 dimensions, we include an apparent horizon in the 

, initial data and prove the Riemannian Penrose inequality in a wide number of 

cases by an elementary argument. In the case of 4 -I- 1 dimensions we obtain 
the analogue of Brill's formula for initial data admitting a generalised form of 
axisymmetry. Including an apparent horizon in the initial data, the Rieman- 
nian Penrose inequality is again proved for a large class of cases. The results 
may have applications in numerical relativity. 



^ ! 1 Introduction 



Some of the most important results in general relativity, are the various positive 
mass theorems for asymptotically flat and asymptotically Anti-de-Sitter spacetimes 
(cf. 0121 El 0). Since the mass in these spacetimes is either constant, or in the 
case of the Bondi mass, non-increasing, these theorems restrict to some degree, the 
possible evolution of spacetimes. Thus, for example since the mass of an asymp- 
totically flat spacetime is bounded below by zero and this can only be achieved by 
flat spacetime, these theorems are sometimes taken to show that flat spacetime is 
stable. However the mass is a single number, determined only by the asymptotic 
structure of spacetime. It can give no detailed information about local behaviour 
of the spacetime geometry. Moreover one knows from the singularity theorems that 
under some circumstances in which the mass at inflnity is positive, the occurrence 
of spacetime singularities is inevitable. 

This situation contrasts sharply with that for classical field theories in flat space- 
time. In that case the mass M may be expressed as an integral 

M = [ Tood^x (1) 

of a local energy density Too which is typically positive semi-definite, vanishing only 
for the trivial field configuration. For example for a real scalar field (j), 
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(V0r+^(f)vn0), 



where is a positive potential function, with global minimum at the origin of 

field space = 0. 

For suitable potential functions V{4>) the constancy of the energy places powerful 
restrictions on how the field (/) may evolve. In particular, using Sobolev inequalities 
it is possible to show that the evolution is smooth for all times. 

Now because of the Equivalence Principle, there is no well defined, generally 
covariant notion, of energy density in general relativity, and a generally covariant 
expression like is not available. However, this does not mean that one cannot 
hope to End an expression analogous to 0) valid in a particular set of coordinate 
systems. Moreover, if the analogue of Tqq is positive definite then one may be able to 
restrict the evolution of the spacctime, and perhaps even prove long term existence. 

Just such a formula was found by Brill [S] for the special case of axisymmetric 
time-symmetric initial data on MP which are also invariant under reversing the 
azimuthal angle^ . Using his formula. Brill was able to give the first mathematically 
rigorous proof of the positive mass theorem in this restricted, but nevertheless 
physically interesting, special case. 

It is widely appreciated that Brill's result may be extended in a straightforward 
way to maximal data. It is possibly less widely appreciated that the angle reversing 
assumption can also be dropped. We will review and generalise Brill's result to 
this so called weakly axisymmetric case in section 2. More interesting physically is 
the possibility of extending his theorem to include an apparent horizon of area A 
say. In this case one expects, because of the cosmic censorship hypothesis |Sj the 
stronger "isoperimetric" or "Riemannian Penrose" inequality to hold 

with equality only for the Schwarzschild initial data. 

In fact ^ has been proven in the general time-symmetric case by Huisken and 
Ilmanen using the inverse mean curvature flow (see also The purpose of this 
paper is to show that Brill's formula may be extended to incorporate an apparent 
horizon and, at least in a large class of cases to be defined precisely below, to 
establish the isoperimetric inequality for black holes. These results are presented in 
section 3. 

In recent years, the global behaviour of higher dimensional spacetimes, espe- 
cially five-dimensional spacetimes, has attracted a great deal of interest because of 
various applications to String and M-theory, and from particle physics phcnome- 
nologists interested in scenarios with small extra dimensions which might permit 
the production of small black holes in accelerators. It is natural, therefore, to en- 
quire about a possible generalisation of Brill's formula to five dimensions, that is to 
four-dimensional initial data sets. We will obtain such a generalisation in section 4, 
subject to the assumption that the four-dimensional initial data admits an action 
by two commuting circle groups, i.e. by the torus group T^. 

In five spacetime dimensions, the isoperimetric inequality for black holes be- 
comes 

,^ St: f AY 3/A\^ 

As yet, there is not a lot of mathematical evidence for this inequality. The 
inverse-mean curvature method is not available because it makes essential use of 

^that is the metric admits an isometric reversible circle action 
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the two-dimensional Gauss-Bonnet theorem. However some support is given by the 
behaviour of coUapsing shells [5] . In section 5 we will use our five-dimensional Brill 
formula to prove the inequality in a wide number of cases. 

Finally, we comment on possible applications of our findings to numerical relativity. 
Three-dimensional axisymmetric initial data sets containing black holes have been 
studied numerically for quite some time now (see e.g. where the superposition 
of Brill- waves and a black hole is studied). Our results concerning four-dimensional 
initial data, in particular formula 19U|I . may be useful for controlling the accuracy 
of numerical codes within the study of higher-dimensional gravity. 



2 Brill's proof of positive mass 

Consider a maximal axisymmetric initial data set (R'^, gab, Kab) for Einstein's equa- 
tions. This means in particular that g and K satisfy the constraint equations for 
trK = 

R = IGtt^ + KabK"^ , (3) 
D'^Kab^O. (4) 

Here D is the covariant derivative with respect to the metric g and is the Tqq 
component of the energy momentum tensor. Axisymmetry means that the metric 
can be written in the form 

g = qAsdx'^dx^ + {d(j) + Aedx^f , (5) 

with QAB a two dimensional metric on the orbit space of the Killing field ^ and 
the functions X and As not depending on the variable (p. If the metric is strongly 
axisymmetric it enjoys an additional mirror symmetry and A has to vanish. Coor- 
dinates can be found such that 

g = e(-2t^+29) ^ ^^2) ^ p2^-2U ^ ^^^^ ^ ^^^^)2 _ 

This choice of coordinates corresponds to finding a harmonic function on the space 
of orbits, i.e. a solution to the equation 

A,,,P = 0. (7) 

The solution p is unique once we specify conditions at infinity and the z-axis, the 
boundary of the orbit space. Furthermore, to obtain a regular axisymmetric metric, 
the function q has to satisfy 

dq 

q = and — = on the z-axis. (8) 
op 

The first condition is necessary to avoid conical singularities on the axis. Further- 
more we impose that at infinity 

1 5 1 

and — (9) 

for some e > 0. On the other hand, the i-term of U contains the ADM-mass of the 

' r 

metric: 

U^-^ and #C/^^ (10) 
r or 

at infinity. 
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Let us compute the scalar-curvature of the metric The dreibein is 
^ e-'-'+idp , e^^e-^+'^dz , = pe''^ {df/) + A^dp + A^dz) . (11) 
The antisymmetric connection coefficients are 

- (12) 
^\ = (i-f/,p)e^-'e3-ipe^-^M^P,.-^.,p)e^ (13) 
= -C/,.e^-9e3 + ipe^'2«(A,,,-A,,p)ei. (14) 

Note that indices are raised and lowered with the Euclidean metric. The curvature 
two-form is obtained from 

i?^ = du'j + uj\ A uj'^j . (15) 
Since the scalar curvature is obtained from 

R = g'^^Rab = 9°''' {R^alb + a2b + a3b) ~ 2i?i212 + 2i?i3i3 + 2i?2323 , (16) 

we only need the A e^-part of i?^2: the A e'^-part of R^^, and the A e'^-part 
of i?^3 to compute the scalar curvature. We obtain 

i?i2i2 = e2^-2''((C/-g),^^ + (C/-g)J-^pV^-4^(Ap^,-A,,p)\(17) 



P J \P 



(18) 



i?2323 = - (C/,.)'e2^-'^ + (e^-''C/,.)_,e^-^ 

+ ip^e^^-^^ {A,. - A.^.f + {U- ,),^ - C/,p) e^^-^ . 

From (|16|l the scalar curvature then is 

(3)^ ^ 2e2^-^^ (([/ - g),^^ + (C/ - ,),,,) - (A,, - A,,) 

which we can write in the form 



(19) 



(20) 



4<^)i?.e(---^) = -Uu+l (V^)Vi A,-ii| + l,^e-^ (A„. - A.^.f . 

(21) 

The Laplacian and gradient are taken with respect to flat coordinates in M.'^. We 
can integrate this expression over M.'^ to prove positive mass. If we take conditions 
(|H1) , ® J into account the integration yields 



M 

16n 



(22) 
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Defining 



Hab = Aa,_ 



A 



B,A ■ 



(23) 



where A, B G {1,2} denote indices on the orbit space with metric qab defined 
above, and noting that X = pe~^ , we can rewrite (I22II in the neat form 



M 



IQtt 



^^^R+ -X'HabH' 



{VUfd^x. (24) 



Since ^^^R = IQi:^ + KabK°-^ > this proves positive mass for the type of initial 
data considered above. 



3 The Generalisation to Black Holes 

In this section we will generalise Brill's formula to initial data leading to black holes. 
To determine the event horizon of a black hole in a four-dimensional spacetime it 
is necessary to know its entire global structure, i.e. there is no direct way to sec an 
event horizon in the initial data. On the other hand, if we consider initial data whose 
maximal development admits a complete null-infinity, an apparent horizon present 
in the initial data will always lead to a regular event horizon in the development (cf. 

Moreover, in this case there will exist on the initial data slice a component 
of the event horizon outside the apparent horizon or coinciding with it. We note 
that for time-symmetric and also for t — 0-symmetric initial-data^ the notion of 
an apparent horizon and a minimal surface coincide. We will restrict to strongly 
axisymmctric data in the following. 

3.1 Static initial data: Weyl's class 

Consider Weyl's family of axisymmctric static vacuum solutions to Einstein's equa- 
tions: 



,2<;(p,z) 



with the functions U and q satisfying the equations 



AC/ 
dpq 



0, 



(25) 



(26) 



d,q = 2pdpUd,U. 



Note that the natural spatial slices of H25|) arc precisely the ones under consideration 
© as initial data slices. If the topology of the initial data slice is trivial, i.e. M"^, 
then by Liouville's theorem, U = (U vanishes at infinity due to asymptotic flat- 
ness). The only static solution in Weyl's class with M.'^ topology is Minkowski space. 
To get non-trivial static initial data, wc will need to allow non-trivial topology. The 
standard example is provided by the Schwarzschild data discussed in the next sub- 
section. Prior to this we explicitly check the validity of Brill's equation 121|) for the 
spatial slices in Weyl's class. Using the equations (|26|) one easily shows that 



0. 



(27) 



Since also AC/ 
holds. 



and i? = for the spatial slices in Weyl's class, equation H21(l 



^At — (/)-symmetric initial data set is defined to be an initial data set leading to a t — 0- 
symmetric spacetime, i.e. a spacetime symmetric under a simultaneous change of sign of the time 
coordinate t and the azimuthal coordinate 4>- 
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3.2 Schwarzschild initial data 

The function U in (|25|l leading to Schwarzschild data can be obtained by calculating 
the Newtonian potential of a rod with unit mass density occupying the z-axis from 
— M to M, as depicted in the figure below. 




We note the relations 

rl = p'^ + {z + MY 
rl = p' + iz-Mf 

The Newtonian potential at the point (p, z) is 



Us{p,z) 



M 



1 



-^log 



^p^ + {z^sf 

M + z + n 



dz 



-M + z + r2 

The function q is then found from the equations H26(l to be 
qs{p,z) = -\og{ 



(28) 
(29) 



(30) 



1, fri+r2+2M 

— log 

2 ^\ri+r2-2M 



(31) 



It will be useful to work in spheroidal coordinates (A, p). They are defined by the 
relations 



with inverse 



and hence 



2A 
2/xM 



- »'l + ?'2 , 

- ri-r2. 
A + pM , 



r-2 = A - p.M , 



p' = {X' - Af-) {I - p') 
z = A/i . 



(32) 
(33) 

(34) 
(35) 



(36) 
(37) 
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Note that e [—1,1]. The transformations 



d d / 1-Ai2 a , d B X^-A'P d 



and hence 



d V(A2 - AP) (1 - A^^) 

^ = p _ (A5a - A^a,) , (39) 

^ = p_),2^j2 (a' - M') d, + X{l- p^) d,) , (40) 
will turn out to be useful later. In spheroidal coordinates we have 

Us^-bo,(^). (41, 



2 " V A- M 



1 / , 

gg= olog x2 ..2. .2 ■ (42) 



One finds that 
dp'^ + dz'^ + p'^dcj)'^ 

and hence the metric (|25|l becomes 



2 ''yX^-p^M^ 



(43) 



- ^^^^ - ^^^^ - ^^^^ - - - ^^^)) • 

To see that this is in fact the Schwarzschild metric, wc finally transform to coordi- 
nates 

A = r - M , (45) 
p = cose, (46) 



to find 



2M\ 2M\ ^ 



ds^ ^[l j di^ ^[l j dr"^ - {dB^ + sin^ Qd4?) . (47) 

Summary: We found that an initial data slice for the Schwarzschild metric can be 
retrieved from the ansatz © for arbitrary axisymmetric initial data by choosing U 
and q as given in H3U|) and H31|l . 

Remark: This is just one possible choice. Another example is provided by choos- 
ing g = and e~^^ = (l -I- ^) , which leads to a conformally fiat, time-symmetric 
Kruskal slice in the Schwarzschild spacetime. However, these spatial slices are not 
the ones used in the Weyl ansatz (|25|l . 

We see that that in the Weyl-form of the metric the Schwarzschild horizon is rep- 
resented by a rod in the initial data. On the rod X = AI and p, runs from —1 to 1. 
The rod has area 

/ e"2^^+«^pdzd0= / / Ul^dpdct) ^ 16ttAI^ (48) 

Jo J -M Jo J -1 
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and the normal derivative 



- VT~^ / / (1 + pdp {-2Us + qs)) dzd<t> = 



rod 



dp 

lULL 

(49) 

vanishes on the rod showing that it is indeed a minimal surface. Equation (|49|l is 
most easily seen in spheroidal coordinates, where one computes (now using (|39() and 



and hence 



as well as 



and therefore 



a 



(51) 



p-^qs 

op 



= 1, (53) 

rod 



from which (|49|l follows. 



3.3 Non-static data: Deformations of Schwarzschild 

Consider now a deformation of the Schwarzschild data. 

U = Us + U and q ^ qs + q , (54) 
where we make the following assumptions about the deformation {U,q): 

1. No new minimal surface should be created by the deformation. (Note however, 
that the deformation does not have to be small in general.) 

2. The rod should remain a minimal surface. This implies that dp (2U -qj =0 
will hold on the rod. 

3. The area of the rod should be unaltered under the deformation: 
A = SttM J^Ij c:ij> (^-~2U + qj dz ^ WttM"^. 

We will eventually need the stronger assumption 

3'. The deformation {U,q) has support only outside the rod, 
i.e. in particular U ^ q = Q on the rod. 

Assumption 3' implies 2 and 3. Near infinity we assume the asymptotics 

-AA/ ,1 
U — and q ~ , (55) 

and that derivatives lower the power of r by one. Integrating H21|) now yields the 
formula 

+ i IIm + h Irod U hn) VUs ^ ^ ?7 A [/s , (56) 
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where the last term vanishes because Us is harmonic. 
If we introduce the quantity 



E 



1 

16^ 



Re 



-2U+2q 



vol 



1 

8^ 



(vuY>o, 



we can write 
AM E 



M 



u 



dz 



p=0 



vol 



M g 
P 



dp 



U 



dz 



(57) 



(58) 



p=0 



describing the change of mass under the deformation considered. The last term 
of H58|) vanishes by assumption 2. Unfortunately, assumption 3 will only give an 
upper bound on the second term of H58|l . since by assumption 3 and the inequality 
> 1 + X we have 



A = IGttA/^ = 8ttM 
and therefore 



M 



M 



exp ( ~2U + q]dz> 16ttM 



8ttM 



M 



M 



M 



U - -q] dz <0. 



-2U + q] dz, 
(59) 

(60) 



Consequently, the second term in (|58|1 could in principle spoil the inequality we 
want to obtain. On the other hand, assuming 3', equation H58|l immediately yields 
a version of the Riemannian Penrose inequality: 

Lemma 3.1: For axisymmetric deformations {U,q) of Schwarzschild initial-data 
satisfying assumptions 1 and 3' the Riemannian-Penrosc inequality 



M > 



A 

16^ 



(61) 



holds. 



Proof. If the deformation has support only outside the rod, both integral terms 
in (|58|l vanish. Consequently, the mass increases by E while the area remains 
unaltered. □ 

One would really like to show more, namely that the inequality H61(l does not 
only hold for arbitrary axisymmetric deformations of Schwarzschild away from the 
horizon but for general axisymmetric initial data containing a single horizon. One 
would like to show the following statement, referred to in the future as "RPI". 



Riemann-Penrose-inequality (RPI): Given any asymptotically- flat, axisymmetric 
initial data with apparent horizon^ of area A. Then (|61|l holds with equality if and 
only if the initial data is Schwarzschild. 

Though we will not be able to establish the above statement by our methods in 
this paper, wc will indicate in the following how one could proceed and what final 
obstacle one encounters. 

To establish RPI it would suffice to show that any axisymmetric initial data can 
be written in Weyl-form as a "deformed Schwarzschild data" with the deformation 
satisfying —2U + q = on the rod. The following lemma shows that it is at least 

''Wc remind the reader that for t — (f> symmetric initial-data the notion of an apparent horizon 
and a minimal surface coincide. 
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always possible to transform the initial data into Wcyl form Ht)2|) using a harmonic 
map. However, we cannot prove that the functions U and q obtained by bringing 
the data into Weyl-form have to agree with their corresponding Schwarzschild val- 
ues on the rod. 



Lemma 3.2: Any axisymmetric initial data set with single connected apparent hori- 
zon of area A, can - in some coordinates - be written in the form 

such that the minimal surface is a rod of area A on the z-axis between —M and M 
for some AI. 

Proof. Given any axisymmetric initial data with apparent horizon of area A, we can 
construct harmonic coordinates (p, z) by solving a particular conformal mapping 
problem. We have to solve 

Ap(p,z)=0 (63) 

with boundary condition p = on the p-axis and on the apparent horizon T and p 
should tend to the distance from the axis, i.e. p, at infinity. 

AP 




This is a classical Dirichlet-problem with Lipschitz boundary, admitting a unique 
solution. The harmonic analogue of p, z, follows from the Cauchy-Riemann equa- 
tions 

dpZ -d^p, (64) 
dzz = dpp. (65) 

The coordinate transformation (p, z) (p, z) maps the minimal surface to a rod 
on the z-axis, lying between —M and M, say. In these coordinates, the metric can 
be written as 

g = e(-2t^+29) (dp2 ^^2) _^ p^e-^'^dcl)^ (66) 

for some functions U and q, unique up to a constant, which is determined by 
asymptotic flatness. □ 

Now given any axisymmetric initial data with apparent horizon of area A, we 
can apply Lemma 3.2, mapping the minimal surface to a rod between —M and Af , 
say, and obtaining unique functions U and q, which are obviously singular on the 
rod. If we split them according to 

U^Us + U and q = qs + q (67) 

where we choose Us, qs as in and ||2U, we find that assumption 2 is satisfied 
and the second term in (|58|l vanishes. Unfortunately we do not know anything 
about the behaviour of U, q (and hence about the first term in l(55)l') on the rod. 
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One way to prove RPI would therefore be to show that the functions U and q 
obtained imphcitly by Lemma 3.2 satisfy the equation 

-2U + q = -2Us + qs (68) 

on the rod. This would mean that the functions U and q, however complicated they 
look like globally, always admit the same blow-up behaviour as in the Schwarzschild 
case, when restricted to the rod. However, we have been unable to show that this 
follows from the harmonic map of Lemma 3.2 and therefore the general case remains 
unproved by our methods. 



4 Generalisation to 4+1 dimensions: Regular data 

In this section we consider regular four-dimensional initial-data for the five-dimensional 
theory of Einstein-gravity. The initial data should admit two Killing vectors with 
t/(l)-orbits. 

4.1 Preliminaries 

The general metric has the form 
g = qABdx^dx^ + Y'^da^ + Z^dad[3 + X'^dfi^ 

+ e^^^''--)^^ sin^ ujda^ + e'^^P^^^p' cos^ udf}^ . 

Without the cross-term W we would have rectangular tori. In general \W\ < 1 
must hold to ensure that the metric on the torus is positive definite. Note that the 
choice f = U = V = W = corresponds to the flat metric 

S = dx'^ + dy^ + du^ + dv^ (70) 

written in the coordinates 

X + iy ^ p sin we'" u + iv — p cos Loe^^ . (71) 
The range of coordinates for (|()^ is 

aG(0,27r] /3e(0,27r] lu e (0,^) . (72) 

There are coordinate singularities at w = and a; = ^ . To avoid conical singularities 
we have to assume certain behaviour of the functions /, U, V, W. We write the metric 
in the orthogonal coordinates given above 

2f+U+Vy2 ^2U^: 




> _|_ ^2 J.2 ^ y2 

■^2f+U+V g2V 







^2f+U+V^2 




V? + 


+ 


1 xvW 





(73) 



x^ + y^)iu^ 
uxWe^+^ 



^(a;2 + y'^)[u'^ + v"^) 
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Regularity (avoidance of conical singularities) requires the cross terms to vanish on 
the axes of the Clifford tori, i.e. 



2f-U + V =0 for u 
2f + U-V =0 for Lu- 



(74) 



and 



W for 



and LU 



TT 

2 



The functions / and W should approach zero for p — s- oo like / 



1 

„2 + t , 



(75) 
and 



W 



with e,S > in order not to contribute to the ADM mass. (Again, 



taking derivatives should lower the order of p by one.) The latter is encoded in the 
lowest order term of U and V: 



2M 



and 



p2 



2M 



(76) 



for p — > oo, where we have defined fi = to keep the notation tidy. 



4.2 Proof of positive mass: Regular data, W=0 

In this subsection we perform the proof of positive mass for the case W = 0, i.e. no 
cross-term in a and /3. Geometrically this corresponds to an additional Z2 symme- 
try. One can think of the Clifford tori as being rectangular in this case. The data 
are supposed to be regular: They have topology and do not contain apparent 
horizons. 

We can perform two Kaluza-Klein reductions of the metric (|69|l , first along ^ , 
then along S^- The result is 



Y 



XY 



(77) 



Note that 



A 



(4) 
flat 



e 



2/-C/-V A (2) 

'^flat ' 

3d 1 92 



r U V 
•^+2 + 2 



1 



d_ 

p dp duP' \ sin oj cos uj J duj 



(78) 
(79) 
(80) 



We compute the terms on the right hand side of (|77|l . We will omit the index "flat" 
from now on since all gradients and laplacians refer to the flat metric. 



X 

A(2)r 

Y 

VXVY 
XY 



A(^V + (Kp)' + 



v.. 



I COSLJ 



smw 



AWl/+([/,)V% + % 
p^ p^ 



sm Lli COS Lli 
COS Lli sin Lo 



sm Lo cos U! 
1 I cos a; 



smu) 



smoj 

COSLJ 



^2f+U+V 



d 



AW -2 A _ 1 

pdp p^ \ sinw coscij / duj 



^+2 + 2 



(81) 
(82) 
(83) 

(84) 
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Hence equation (|77|l becomes 



.1(4)^ . e2/+c^+v = ^AU+-AV + Af 
2 2 2 



1 



2p2 sin a; cos w 



{u - V) 



(85) 



1 



smw 
cos a; 



We now integrate over M^. The term on the left will be manifestly negative since 
the Ricci-scalar is positive by the constraint equations for maximal data. The first 
two terms on the right hand side will give —SMtt by the assumption H76|l . The 
term A/ vanishes upon integration because dpf decays like -^jtjt- The terms in the 
second line are manifestly positive. Finally, we will now show that the terms in the 
last two lines cancel when integrated over M.'^. We have 



47r2.1 y" pd^ [U - V) dpdco - / p 

and 



{U-V) (^p,^)-{U-V){p,0)]dp (86) 



-47r^ • 2 / p^ f^p sin uj cos ujdpduj = — 47r^ • 2 



p—oo 

sin uj cos toduj 
P=o (87) 



Att ■ 4: pf sin UJ cos uidpdcu , 



and 



47r^ J p (sin^ u) — cos^ oj) d^jdpdio = 
/ 1^/9, — ^ + / (p, 0) dp — Att'^ ■ 4 / pf sin UJ cos ujdpduj . 



An' 



The last term of (|88f) cancels with the last term of (|87f) . The sum of the expression 
(I86|l and the first term of H88() vanishes due to the regularity conditions (|74|l imposed 
on the axes. Finally, the first term of H87|l vanishes because / falls of faster than ^ 
at infinity as assumed above. 
We can summarise the result in the formula 



M=— ( Wi?.e2/+c^+Vd4x+ 
167r 7r4 



1 

8^ 



i(vc/ + vy)' 



(89) 



which clearly shows M > 0. 
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4.3 The case 7^ 



In the W =/= case the calculations are much more involved but straightforward. 
The analog of (|85|1 reads 



_i(4)^ . e2f+u+v = ^AU + -AV + Af 
2 2 2 

+ I + i (S/Vf + i (V(7 + vy)' 

+ TT^-^ {U - V) 

Zp^ smuj cos a; 

2 If sinuj coauj\ 

-dpf + : d^f 

" \ coscj suicj / 




3 {WWf _ W^jVWf _ W IdW 

4 1-1^2 ~ (1 _ ^2)2 " 1-M^2--^ 



U,^ W,^ cos(2w) , 1 1 



2/9^ 1 — W"^ \ sin w cos u sin a; cos w W sin w cos w 2 sin^ uj cos^ , 

(90) 

We already know how to deal with the terms of the first four lines from the previous 
section. We will now use partial integration to handle the terms in the fifth line 
including the function W only. We observe that 

W AW ilVWf W^iVWf W ldW\ o . , , , 

— - — p sm Lo cos Lodpdujdadp 



1-1^2 4 1 - (2 _ VK2)2 l-W"^ p dp 

'l{VWf W^{VWf , 2^ 1 



4 1 - (1 _ 1^2)2 o V / p 



— log (1 — W^) ) p^ sin UJ cos Lodpduj dad P . 



(91) 



All terms on the right are manifestly non-negative. We will need the first and 
the second term of (I91II in the following to make other expressions manifestly non- 
negative. Next we collect the terms of (I90|) involving a derivative oi U oi V in p, 
borrowing also one of the manifestly non-negative terms we just obtained in H91|l . 
We will show non-negativity of the following expression: 



W 1 

iUpW^, + VM + --—2 {Up - Vpf > . 



(92) 



2 1-1^2 y--p--p ■,P',p^ Al~W^ 



(The last term in the first line of (|92(l comes from (|91|l .') Defining a = U^p, b — V 



c = ^-w^ ^® need to show 



3 3 

+ b'^ + + ab- -be - -ac > . (93) 

For this we simply write the left hand side as 

i((a-6)2+c2) + J(a + 6-c)^>0. (94) 

Next we turn to the terms of (|90() involving an uj derivative and the term just 
involving trigonometric functions. Collecting these terms (again borrowing two 
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manifestly negative terms of (|91|l ') we have to show non-negativity of the expression 



1 

V 1 - 



C/,^ - - ^ cos(2c^) + i 

W sm uo cos ^ 

1 (VF,^)'cos2(2a.) , 1 



4p2 



1 - 



4p2 1 - Vl^2 



1 W^^VF,,, 



(95) 



cos(2w) 



p2(l- VF2) 



p2(l-M^2)2 

(c/,..w^,^ + ^^w^,^) >o. 



(Note that the second term in the second hne and the last term in the third line of 
(I95|l are taken from (|91|l .) The terms in the first and second line of the expression 
(I95|l are already manifestly non-negative. For the remaining terms we define A — 
B = and C = and show that 



1 



1 



A'' + B^ + + AB - AC y-- - cos(2w) j -BC\- + - cos(2tj) ) > . (96) 
To see this, we just write the left hand side as 

3 . . ^ "^2 1 „ , ^ ^^2 1 



^{A + B-Cf + ^{A-B + C cos{2uj)f + ^C^ (l - cos2(2tj)) > , 
making it manifestly non-negative. We finally obtain the formula 



(97) 



M=— W^.g2/+c/+y^4^ 



1 

167r 
1 

8^ 



iog(i -w^)^ + + 



+ \{{a^hf + c^) + \{a + b-cf + \{A + B~Cf 



1 



1 



+ - 



-{A-B + Ccos{2uj)Y + -C^ (1 - cos^{2uj)) 

1 1^2 / I 

[A-B--^ cos(2cj) + 

\ Wp psinwcostj 

(a - b f ] d^x . 



41 
1 



41 -W^ 



(98) 



5 Generalisation to 4+1 dimensions: Black holes 

In this section we generalise the Brill formula to include black holes and derive a 
Riemann-Penrose inequality for the initial data (cf . Lemma 5.1). We restrict to 
the case of rectangular tori. 



5.1 The generalised Weyl-Class 

The analog of the four-dimensional Weyl-class for orthogonal tori was obtained in 
|12| . Here we will use a slightly different representation and state the relevant 
formulae in both toroidal and standard r — z-orbit space coordinates. The metric 
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g = _^-2U{p,u.)-2VM^^l ^ ^2f(p,u.) + U(p,u.)+V(p,u) ^^p2 ^ p2^^2^ ^ 



(99) 



with Einstein equations 

^flat ^ - 



3 I 1 0^ 
p dp ~p^ diu'-^ 



1 / cos LJ _ sin uj '\ 
p-- V sin uj cos a; / duj 



U{r,iu)^0, (100) 



and 



--dpf 
P 



4^ + -^ + ^w^ + ^(^-—)w-]uir,u;)^0, (101) 

ap^ P f^P P ouj^ p^ \ sin UJ cos cj / au; v ' y ' v / 

4 cos(2a;) 



sin(2a;) 



{dpuy + dpUdpV + {dpvy 

{dujUf + d^Ud^V + {dpVf 
^ (dujU-d^V), 



costj smtj 



(102) 



p sin(2w) 



- {2dpUd^U + 2dpVd^V + dpVd^U + dpUd^V) 



1 



psin(2a;) 



[dpV - dpU) . 



(103) 



The relation of the toroidal variables p and uj to the the standard r-z orbit-space 
variables used in |12| is given by 



r — —p^ sin (2a;) 
z = —p^ cos (2cj) 
such that in z — r-coordinates the metric (|99ll reads 

g ^ _^~2Uir,z)-2V{r,z)j^^2 ^ ^2f {r,z) + U {r,z) + V {r,z)-\oe{2V^F^l^) ^^^2 ^ ^^2-j 

_^g2C/(r,^) ^ v/r2 + da^ + e2^('^^) (^z + y^r^ + z^^ dp^ . 
The equations for the functions U and V translate to 



2^7-2 + z2 



2\/r^+z^ 



dl + -a, + dl 

T 



di + -a, + di 

r 



U{r, z) = 
V{r,z) = 



(104) 
(105) 

(106) 

(107) 
(108) 



involving the three-dimensional flat Laplacian just like in the three-dimensional 
case. It is clear how to translate the equations (|l()2(l and (| 10311 into — z-coordinates. 

Using equations H100|l . (|101() . p02|l . H103() and the fact that the Ricci-scalar of 
the initial data slice vanishes because of the time-symmetry of the slice and vacuum, 
one easily checks that the differential Brill-formula (|85|l holds for the spatial slices 
in the five-dimensional Weyl-class H99|l . 
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5.2 Schwarzschild-Tangherlini initial data 

A conformally flat timc-symmctric initial data slice for the Schwarzschild-Tangherlini 
[T5| metric can be obtained by choosing f = W = and 

^\ , / 2M 



i; = f = log I 1 + ^ = log 1 + — , ) (109) 



in the ansatz H69|l . This would lead to a spatial slice of the form 

5 = + {dp^ + P^diol) , (110) 

for which we find a minimal surface at p ^ ^i^h area 

A^i&^{^-fy . (Ill) 

In particular, this slice covers both asymptotically flat ends. 

However, this slice is not part of the slicing implicit in the Weyl-ansatz (|99|l . The 
Schwarzschild metric in Weyl-coordinates is given by a certain rod-representation, 
just as in the four-dimensional case. We will now turn to this representation, giving 
all formulae in r — z coordinates, i.e. the metric functions for the ansatz (|106|l . 
They can easily be translated to toroidal coordinates using (|104() . The harmonic 
functions U and V are 



^ — z + V + z^ 



1 /i + z + \/(m + z)^ -f 

^s--^^og . , . (113) 

^ z 4- V f + z^ 

The function Us corresponds to the Newtonian potential induced by a rod lying 
on the z-axis in the interval [0,/i] whereas Vs corresponds to a rod in [— //,0]. 
Asymptotically, 

Us ^ , ^ = ii Vs'- . ^ = ^ (114) 
as it should be. The function fs is found to be given by the expression 



/. = -^-^4log(V^)+ilog(- 
where we have defined the distances from the rod 



2 



P 



rir2 



(115) 



Ti = \j {z + pf + r"^ r2 ^ \J {z - p)^ + r"^ (116) 

analogously to the four-dimensional case. The function fs is unique up to a con- 
stant, which is determined by asymptotic flatness. Asymptotically, 

/s-2(r2 + z2)-7r' (117) 
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as required. Inserting Us, Vs and fs into (|106|l will yield the Sehwarzschild met- 
ric. Transforming to toroidal coordinates (|104|l and then applying the coordinate 
transformation (cf. \V2\ ) 



w = - arctan 



1 



{R* - + cos (261)) 5 



iii 

tan (20) 



(118) 



(119) 



will finally lead to the Sehwarzschild metric in the standard {t, R,9,a, (3) coordi- 
nates: 



9S 



1 - ^\ dR' + R' 



f + sin^ eda^ + cos^ 9dp^) 



where /i = as defined above. 



(120) 



3tv 



As a check we compute the area of the rod using (|99|) . The right part of the rod, 
located on the z-axis between and fi, has area 



A 



Js + \Us + %Vs J. 



p cos Lo sin uj \ dpdadp 



(121) 



Note that this will only give half of the area of the full rod, since we also need to 
take the left part (located at w = ^) into account. Inserting the functions Us, Vs 
and fs defined above, we find that 



and hence 



gfs + lUs+lVs p2 coswsinw 



(122) 



(123) 



which - if we take both halves of the rod into account - is consistent with (|lll|l . 
the area obtained for the conformally fiat slice. 

Finally, using H121|l one checks that the rod between —p and p is indeed a minimal 
surface, namely that 

l^A = ±p-^A = (124) 
p OLo or 

holds on the rod. 

5.3 Brill- formula 

Just as in the four-dimensional case, we can consider deformations of the Sehwarzschild 
initial data: 



U^Us + U V = Vs + V f^fs + f, 
where we assume the following asymptotics near infinity 



U 



2AM 
Snp^ 



V 



2AM 



/ 



1 



(125) 



(126) 
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and that taking derivatives lowers the order of p by 1. Integrating the differential 
Brill-formula H85|) will now lead to 



Idtt / Stt 



1 

r 



— TT 

2 



3 ~ 



psinwcoswSij -J7 + -y + / \Pi7:]dp 



/27Z 



/3 ~ 3 ~ 
psinwcos^a^ ( -[/+ -y + / ) {p,0)dp 



(127) 



sm(2to') 



2 = cos(2(jj) 




which can be translated to ?' — z coordinates, where we find the expression familiar 
from the three-dimensional case (cf. equation (|58|l 'l 



AM = — / (4)^g2/+c/+y ^ ^ 



dr \2 



r=0 

dz . 



(128) 



To obtain these formulae one uses 

-n- V 



rod 



ld_ 

p dui 



dr ' 



(129) 



where the upper (lower) sign stands for right (left) part of the rod, = (cj = f ) 
The vector —n is the unit-normal vector on the rod. Furthermore, 



4us 
or 



— 1 on the right part of the rod 
else 
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d 
or 



— 1 on the left part of the rod 
else 



From (|127|l one infers the following analogue of the Lemma 3. 1 : 

Lemma 5. 1 : For T^-symmetric deformations of Schwarzschild-Tangherlini initial 
data (/, U, V) admitting support only outside the minimal surface, i.e. U = V = 
/ = on the rod, the mass of the new initial data will increase, whilst the area will 
stay the same. Hence 

holds for deformations away from the minimal surface. 

To prove the four-dimensional Riemannian-Penrose inequality in full generality in 
the T^-symmetric setting it remains to show that any biaxisymmetric initial data 
can be written in deformed Schwarzschild form such that |C/ + + f ~ always 
holds on the rod, leading us to the same obstacles encountered in section 3. 

6 Conclusion 

We have presented various generalisations of an argument used by Brill to prove 
positive mass of three-dimensional, axisymmetric, time-symmetric, regular initial- 
data for Einstein's equations. We first noted that the symmetry assumptions can 
be slightly weakened. Furthermore an apparent horizon was introduced in the data, 
which was seen to be represented by a rod, when the data was put into Weyl's 
form. A simple calculation then established a version of the Riemannian Penrose 
inequality, namely that all axisymmetric deformations (not necessarily small) of 
the Schwarzschild metric away from the horizon will increase the mass while leav- 
ing the area invariant. We indicated how we hope to proceed with an elementary 
proof of the general case in this axisymmetric class. While the Penrose-inequality 
has been proven for the time-symmetric case in great generality, there is so far not 
much mathematical evidence that it will hold in higher dimensions. Here we were 
able to show the version of the Riemannian Penrose inequality mentioned above for 
four-dimensional initial data admitting an action by the torus group T^. For the 
proof we derived the generalisation of Brill's formula for 4+1 dimensions and also 
obtained the positive mass theorem for data with trivial topology. 

Although the primary objective will be the removal of the gap preventing us to 
prove the Riemannian Penrose inequality in full generality using our method of 
rods, the techniques used in this paper should also allow generalisation to further 
interesting cases, which we hope to address in the future. On the one hand, one 
could include more than one apparent horizon and try to extract quantitative infor- 
mation about the interaction energies and the gravitational radiation produced in 
the development of the data (see m| for an attempt in this direction) . On the other 
hand, one could enquire about the inclusion of a cosmological constant or a study 
of the electrovacuum case. In the latter case, one may hope to prove the stronger 
inequality 

as suggested in . All generalisations mentioned can be attempted for both 3-1-1 
and 4+1 dimensions. 
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